Hot Electrons in Germanium and Ohm's Law 

By W. SHOCKLEY 

The data of E. J. Ryder on the mobility of electrons in electric fields up to 
40,000 volts per cm are analyzed. The mobility decreases many fold due to the 
influence of scattering by optical modes and due to increases of electron energy. 
It is estimated that electron "temperatures" as high as 4000°K have been pro- 
duced in specimens having temperatures of atomic vibration of 300° K. The 
critical drift velocity above which there are deviations from Ohm's law is about 
2.6 X 10* cm/sec. This is three times higher than the elementary theory 
and an explanation in terms of complex energy surfaces is proposed. 
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1. Introduction: Fundamental Deviations from Ohm's Law 

THE starting point of many branches of physics is a linear relation. 
Among the most prominent of these are Hooke's law, which relates 
stress and strain for solid bodies, Newton's second law of motion F = 
ma and Ohm's law. In all of these cases, the linear relation is only an ap- 
proximation that may be regarded as the first term in a Taylor's expansion 
of the functional relationship between the two variables. Important physi- 
cal principles are brought to attention when the nonlinear range is reached. 
Of the three laws mentioned, Newton's is, of course, the one in which the 
failure of linearity is the most significant representing as it does the en- 
trance of relativistic effects into the laws of motion. 

The failure of Hooke's law may be of either a primary or secondary form. 
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If a solid contains voids, then under a certain pressure it will crumble and 
fill the voids. This is a secondary effect. If the sample is homogeneous, 
however, high pressures will produce fundamental deviations from Hooke's 
law, these deviations arising from the nonlinearity of the forces between 
atoms. Studies of these nonlinear effects by Bridgman have, among other 
things, put on a firmer basis the understanding of the forces between ions 
in ionic crystals and the pressures of electron gases in metals. 

Deviations from Ohm's law for electronic conduction in semiconductors 
are almost the rule rather than the exception, but the most familiar cases 
are secondary rather than primary. The primary linear relation for the 
conduction process is that between the drift velocity of an electron, or hole, 
and the electric field that drives it. This relationship is 

v d = M „ (r) e, (i.i) 

where the mobility no(T) is a function of the temperature T of the specimen. 
By a fundamental deviation from Ohm's law we shall mean a deviation in this 
linear relationship arising from the largeness of E rather than other causes. 

Thermistor action is typical of a secondary deviation from Ohm's law. 
A thermistor is usually a two-terminal circuit element in which the current 
flows through an electronic semiconductor. The semiconductor has the 
property that its resistance decreases rapidly as the temperature increases; 
and the physical basis for this decrease is an increase in the number of con- 
ducting electrons (or holes or both) with increasing temperature. The 
passage of current heats the thermistor and its resistance changes; conse- 
quently the linear relation between current and voltage fails and in fact 
there may result a decrease of voltage with increasing current so that a 
differential negative resistance is observed. The electric fields are so low, 
however, that equation (1.1) is valid provided the dependence of mo on the 
temperature is taken into account. An experimental proof that no funda- 
mental deviation of Ohm's law occurs is furnished by applying a small 
a-c. test signal on top of a d-c. bias that produces heating. If the frequency 
is much higher than the thermal relaxation rate, the a-c. resistance is found 
to be simply that expected for the observed temperature. 

The principal nonlinearities of crystal rectifiers, or varistors, and of tran- 
sistors are also secondary and are associated with changing numbers of 
current carriers. 

In this article we shall discuss some experimental evidence of funda- 
mental deviations in Ohm's law for electrons in w-type germanium obtained 
by E. J. Ryder of Bell Telephone Laboratories. 1 We shall describe his ex- 

l E. J. Ryder and W. Shockley, Phys. Rev. 81, 139 (1951). 
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perimental techniques and results briefly in the next section and shall then 
present some aspects of the quantitative theory that explains them, leaving 
the bulk of the mathematical manipulations for the appendices. 

Before discussing Ryder's results, we may indicate why his procedure 
succeeded whereas previous attempts, of which there have apparently been 
a number, largely unpublished, have failed. Ryder's work takes advantage 
of three factors: (1) the availability of electrical pulses of microsecond dura- 
tion, (2) the high resistivity of germanium, and (3) the high mobility of 
electrons in germanium. Because of (3), it is possible to deliver energy to 
electrons at relatively high rates by electric fields. In effect this "heats" 
the electrons above the temperature of the crystal and lowers their mobility. 
The generalized equation is then 

v d = p(T, E)E (1.2) 

where the fact that m depends on E represents the fundamental nonlinear- 
ity. We shall show that Ryder's techniques raise the "temperature" of the 
electrons by a factor of about thirteen fold to above 4000°K. Since the re- 
sistivity is high, say 10 ohm cm, the power delivered to the specimen is 
sufficiently low that the heating in one pulse is negligible. The pulse repeti- 
tion rate is then kept so low that accumulated heat is negligible also. 

These conditions are enormously more favorable than those met with 
in metals. In a metal the average electron energy is several electron volts; 
in order to double this energy, each electron would acquire an added energy 
roughly equal to the cohesive energy per atom of the crystal. Furthermore, 
in a metal there is about one conduction electron per atom, compared with 
10 -7 per atom in Ryder's samples. Thus the stored energy due to "hot" 
electrons in a metal would be enough to vaporize it, whereas in germanium, 
or a similar semiconductor, a temperature of 10,000°K for the electrons 
would be enough to raise the crystal less than 0.01°K. From this reasoning 
it appears that it will be extremely difficult, if not impossible, to produce 
significant fundamental deviations from Ohm's law in metals and certainly 
impossible to produce effects of the magnitude described below. 

It should be pointed out that the behavior of electrons in crystals in 
fields so high that equation (1) fails have been subject to both experimental 
and theoretical investigation in connection with, dielectric breakdown. 2 
The work does not apply to cases in which the specimens obey Ohm's law 
at low fields, however, and the experiments do not permit accurate deter- 

2 See, for example, H. Frohlich and F. Seitz, Pliys. Rev. 79, 526 (1950) and F. Seitz' 
Phys. Rev. 76, 1376 (1950). Much of the treatment presented in the Appendices is essen- 
tially equivalent to that given in Seitz. In our Appendices, however, we give much more 
emphasis to the low field case. The Seitz paper also contains a review of the literature 
to which the reader is referred. 



HOT ELECTRONS IN GERMANIUM AND OHM'S LAW 993 

minations of Vj as a function of E. From the theoretical side also the empha- 
sis has been on fields so high that the linear range is neglected so that the 
transition from linear to nonlinear is not stressed. 

The current theories of dielectric breakdown are based on the principle 
of "secondary generation" or "electron multiplication." Thus if an electron 
acquires enough energy from the electric field, it will be capable of pro- 
ducing secondaries by collision with bound electrons, and the repetition of 
this process will lead to an avalanche. Our theory indicates that in ger- 
manium, even at fields as high as 200,000 volts/cm, few electrons will have 
enough energy to produce secondaries. At about those fields, however, 
another phenomenon occurs. 

In 1934 C. Zener 3 proposed that dielectric breakdown was due to a pri- 
mary effect: the field induced generation of hole-electron pairs. His mathe- 
matical theory is similar to that for field emission from cold metal points 
and to that for radioactive decay. It involves the "tunnelling" of electrons 
through regions in which their wave functions are attenuated, rather than 
running, waves. 

Zener's theory does not seem to apply to breakdown; however, it does 
apply to the high electric fields produced in rectifying p-n junctions in ger- 
manium when these are biased in the reverse direction. Under these condi- 
tions fields of the order of 200,000 volts/cm are produced. The mobilities 
of electrons, or holes, in these fields have not been measured. It has been 
shown, 4 however, that secondary production is very small. At these fields 
a sort of "breakdown" effect occurs and above a critical value of the volt- 
age a very rapid increase in current is observed. This current appears to 
be of the nature predicted by Zener. It is stable at a given voltage, has a 
small temperature coefficient and will probably be useful in semiconductor 
analogues of "voltage regulator tubes" and protective devices. 

As is shown in the treatment given in qualitative terms in Section 3 and 
in more detail in the Appendices, the explanation of the fundamental devia- 
tions from Ohm's law is based on the theory of electron waves. The investi- 
gations described in this paper may thus be regarded as furnishing evidence 
for the wave nature of conduction electrons in germanium and are thus re- 
lated to the researches of C. J. Davisson, to whom this volume is dedicated, 
and his collaborator, L. H. Germer. The Davisson-Germer experiments were 
concerned chiefly with electron waves in free space and with high energy 
electrons in crystals. Both of these cases are simpler than that dealt with in 
this paper. Electrons in the conduction band in germanium appear to behave 
as though they were in a multiply refracting medium in which they may have 

3 C. Zener, Proc. Roy. Soc. 145, 523 (1934). 

* K. C. McAfee, E. J. Ryder, VV. Shockley and M. Sparks, Phys. Rev. 83, 650 (1951). 



994 



THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1951 



several velocities of propagation for any specified direction of propagation 
and frequency. It is to be hoped that more detailed analyses of the data ob- 
tained by Ryder, together with quantitative interpretations of certain ob- 
servations of magnetoresistance, may lead to a unique evaluation of the 
"refractive constants" of the medium for the electron waves; this possibility 
is discussed briefly in Section 5. The phenomena in the Zener current range 
afford another new opportunity to study electron waves in crystals. The 
waves involved in this effect are those with energies in the energy gap be- 
tween the conduction band and the valence band; waves in this range have 
received little attention from either the experimental or theoretical side. 

2. E. J. Ryder's Results 

One of germanium's most noted attributes is its ability to give am- 
plification of electrical signals when made into a transistor. The basic phe- 
nomenon for many types of transistors is that of "carrier injection." As is 




Fig. 1-The principles of E. J. Ryder's technique for observing conductivity in high 
electric fields. 

well known, germanium may carry current either by the electron mecha- 
nism in which case it is called w-type germanium, or by the mechanism of 
hole conduction in which case it is called p-type. If a suitably prepared 
electrode is placed on an «-type specimen and current is caused to flow in 
the sense that removes electrons from the specimen, then the process may 
cause "hole injection." In this case in addition to removing conduction 
electrons from the germanium, electrons are removed from the valence 
bonds so that holes are injected. This leads to nonlinear effects because, 
as the current passes through the specimen, the number of carriers in the 
specimen changes and so does its resistivity. 

In order to avoid the secondary deviations from Ohm's law due to carrier 
injection, Ryder has designed specimens of the form shown in Fig. 1. These 
specimens have large ends to which the metal electrodes are attached. The 
resistance arises chiefly from a thin section of the material connecting the 
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large ends. Since the fields at the large ends are small, carrier injection is 
largely suppressed; furthermore, the electric fields are applied for such a 
short time during the pulse that, even if holes were injected at one of the 
ends, they would not have time to reach the narrow section of the bridge 
and modulate its conductivity during the period of the pulse. 

Further causes of non-linearity can arise from inhomogeneities in the 
germanium material itself. For example grain boundaries in polycrystalline 
germanium are known to have added electrical resistance. Difficulties due 
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Fig. 2— Currents and estimated drift velocities deduced from E. J. Ryder's pulse data 
on a specimen of n-type germanium of 2.7 ohm-cm resistivity. [The fact that the numeri- 
cal values of current density and drift velocity have the same digits is a consequence of 
the accident that a = 1/2.7 = 0.37 is almost exactly 10 -4 times the mobility.] 



to inhomogeneity have largely been eliminated in these experiments by 
the use of highly homogeneous single-crystal germanium material fur- 
nished by G. K. Teal and his collaborators. 

Other experimental precautions are necessary, such as assuring a smooth 
polished surface on the filament; if this is not done, apparently holes are 
injected from the surface irregularities of the thin section between the 
large ends. It is also necessary to make corrections for end effects since 
some of the resistance arises within the large blocks themselves. 

Some of the data obtained by Ryder are shown in Fig. 2. The drift ve- 
locity, plotted as ordinate, is not measured directly but is inferred from 
the measured currents through the specimen by the following reasoning: 
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In a specimen at room temperature the drift velocity of electrons is given 
by the equation 

v d = 3600 E cm/sec (2.1) 

when E is expressed in volts per cm. 5 At 100 volts per cm, for example, 
(which is below the non-linear range) the velocity of electrons should be 
3.6 X 10 5 . This establishes the drift velocity scale for room temperature. 
From other measurements of the germanium specimen of Fig. 2, it is con- 
cluded that the number of electrons available for conduction is substan- 
tially independent of temperature down to liquid air temperatures. Conse- 
quently, for this temperature range the drift velocity should be directly 
proportional to the current in the specimen. The other two sets of data 
are accordingly simply scaled in proportion to their currents. 

On the figure we also show extrapolated lines at 45 degrees corresponding 
to Ohm's law. From these we see that decreases in mobility of tenfold or 
more have been produced in these experiments for high field conditions. 

The data for each temperature fall approximately on three lines: The low 
field or Ohm's law region, an intermediate region over which v d is propor- 
tional to .E" 2 and /x is proportional to £~ 1/2 , and a saturation region. The 
break at low fields comes at drift velocity of about 3 X 10 6 cm/sec for all 
three cases. This break, according to theory, should come when the drift 
velocity is several times the speed of sound in germanium, the speed of 
sound being about 5.4 X 10 s cm/sec. The limiting drift velocity at higher 
fields is associated with the energy required to excite a particular type of 
atomic vibration, called an "optical mode." It comes at approximately 
the value of drift velocity predicted by theory. 6 The theoretical curves, 
computed in the appendices, do not break into the sharp line section sug- 
gested on Fig. 2. However, they show the distinct influences of separate 
causes and fit the data reasonably well, as we shall show below in connection 
with Fig. 5. 

3. Theory of Deviations from Ohm's Law 

3a. Electrons in n-Type Germanium 7 

The specimens we shall consider are of «-type germanium and have 
resistivities of several ohm cm. The conductivity arises from the presence 

« J. R. Havnes and W. Shockley, Phys. Rev. 81, 835 (1951). 

6 The observation and explanation of these general features was presented in our 
first publications: E. J. Ryder and W. Shockley Phys. Rev. 81, 139 (1951) and 82, 330 

7 The material under this heading is treated in more detail in the author's book, 
"Electrons and Holes in Semiconductors," D. van Nostrand (1950), Chapter I. This 
book will be referred to subsequently as E and U in S. 
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of donors: chemical impurities such as arsenic or antimony. These donors 
substitute themselves for germanium atoms in the crystal structure, form 
electron-pair bonds with their four neighbors and release their fifth valence 
electron to the conduction band. The density of donors is about 10 l5 /cm 
or one per cube 10 -5 cm = 1000 A on an edge. The donors are fixed positive 
charges and do not move in electric fields. Their charges neutralize those 
of the electrons. The electrostatic energy of interaction between electrons 
and donors leads to a deflection of the electron's motion. For the tempera- 
tures of Fig. 2, however, this effect is unimportant compared to the effect 
of thermal vibrations of the atoms. 

The electrons in the conduction band move in accordance with a wave 
equation. They may, however, be thought of as particles. The justification 
is that, under many experimental conditions, the wave functions will actu- 
ally be wave packets. These wave packets, it can be shown, behave much 
as particles and can be dealt with as particles, at least provided the phe- 
nomena considered do not involve distances smaller than the size of the 
wave packet. 

Under conditions of thermal equilibrium we may think of the 10 1S elec- 
trons in each cubic centimeter as an electron gas with the electrons (as 
wave packets) moving at random with an average kinetic energy of motion 
of (3/2)kT. 

If the atoms of the crystal were held rigidly at rest in a perfectly regular 
crystal structure, an electron wave, and a wave packet too, would be trans- 
mitted through it with no scattering. At 300°K the vibrations are such that 
the wave packet moves for only 2500 A before being scattered. At low tem- 
peratures, the mean free path is longer and at liquid hydrogen tempera- 
tures it is so long that thermal vibrations are less important than the 
fields of the ionized donors. As stated above, however, we may neglect this 
scattering by ions over the temperature range of Fig. 2. 

Before proceeding with the discussion of the interactions of electrons 
and thermal vibration we shall point out that two problems must be solved 
before the dependence of mobility upon electric field can be explained: 

First, the mechanisms of the individual processes must be analyzed. 
This is the basic physical problem. In order to solve it we must apply quan- 
tum mechanics to the model representing the electron moving in the crystal 
and determine the probabilities of various types of transitions and some 
appropriate averages. 

Second, the statistical consequences must be worked out. On the basis 
of the individual processes, the statistics of the assemblage of electrons 
must be analyzed and a steady state solution found. 

The first problem poses the more physical problems and is given the most 
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attention. The second problem is more difficult mathematically. It is given 
only an approximate treatment which is adequate, however, to indicate 
that the solution to the first problem contains the necessary features to 
explain the experiments. 

3b. The Phonons 

We must next consider how to describe the thermal vibrations and to 
evaluate their interactions with the electrons. We shall present only the 
principal results of the mathematical analysis here, leaving the details for 
the appendices. The earliest treatment of thermal vibration in a crystal 
was that of Einstein, who considered each atom to be a separate harmonic 
oscillator. This model was improved on by Debye who treated the crystal 
as an elastic continuum that could support running waves. Debye's method 
is regarded as essentially correct and we, therefore, resolve the atomic 
motions into a set of running waves, or normal modes. There are three 
times as many independent normal modes as there are atoms in the crystal, 
or one per degree of freedom, and any possible atomic motion of the crystal 
may be made up as a sort of Fourier series in these normal modes. 

Each normal mode must be treated as a Planck oscillator and has a 
system of energy levels with values 

(» + l/2)hv (3.1) 

where v is its frequency of vibration. Each quantum of energy is referred 
to as a phonon; if a normal mode makes a transition with 8n = +1, we say 
a phonon has been emitted and if Sn = - 1, we say one has been absorbed. 

The description of the crystal in terms of phonons is in close analogy 
with the description of electromagnetic waves in a cavity in terms of pho- 
tons. For the case of light, the electromagnetic state of the cavity is deter- 
mined by finding the normal modes, which are treated as quantized oscil- 
lators, and transitions with bn = ±1 correspond to photon emission and 
absorption. 

The normal modes for the crystal are unlike those for light. For low 
frequencies the waves are essentially the microscopic transverse and longi- 
tudinal waves of a solid. As the wave length becomes shorter, however, 
the sound velocity varies and there is a limiting minimum wave length 
which is about twice the spacing between atoms. In order to understand 
the energy losses of electrons in high fields, we must consider the role of 
this minimum wave length. For this purpose we shall describe the depend- 
ence of frequency upon wave length for a longitudinal mode. 

Accordingly we consider the frequency of the normal modes correspond- 
ing to a longitudinal wave propagating along a cube axis. Rather than 
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using the wave length X as a variable, we use the wave number or (1/X). 
For long waves the frequency is simply 



„ = c /\ = c(l/X). 



(3.2) 



This corresponds to the straight line portion for low frequencies in Fig. 3' 
This portion extends to a wave length equal to twice the lattice constant a 
of the crystal. 

Figure 3 shows another curve which has a high frequency even for 
(1/X) = or infinite wave length. The presence of this branch of the "vi- 
brational spectrum" is due to the fact that the diamond structure has two 




WAVE NUMBER 



1/2 a 



Fig. 3 — Frequency of longitudinal vibrations in [100] direction in the diamond struc- 
ture. (In this particular direction of propagation the acoustical and optical branches 
join smoothly at the same frequency; for other directions, there is a discontinuity in fre- 
quency. The dependence of v upon 1/X is approximated by a sine wave.) 

atoms per unit cell. (The diamond structure is made of two face centered 
cubic arrays of atoms, juxtaposed so that each atom of one array is cen- 
trally situated in respect to a tetrahedron of four atoms of the other array, 
with which it forms four electron-pair bonds. The unit cell contains one 
atom of each array.) As a consequence of this it is possible to have a vibra- 
tion in which one atom vibrates in the plus x direction while the other atom 
vibrates oppositely and to have this same motion occur in phase in every 
unit cell. Such a vibration is considered to have infinite wave length, since 
every unit cell does the same thing at the same time. It has the highest 
possible frequency since the pattern of motion involves directly opposed 
motions of nearest neighbors. If the motion is modified so as to have dif- 
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ferent phases in adjoining unit cells, and thus to correspond to a finite wave 
length, the frequency drops. 

The opposed motions are referred to as "optical modes" by analogy 
with polar crystals. In a crystal of sodium chloride there is one Cl~ and one 
Na + per unit cell. In the opposed type of motion, ions of like sign move 
one way and opposite to those of the other sign. This relative motion of 
charge polarizes the crystal and phonons of this type of vibration can ab- 
sorb or emit light. Because of this optical activity the mode is termed 
optical. 

The name "optical" is carried over to valence crystals to describe the 
opposed form of motion, although no polarization accompanies the dis- 
placement in the latter case. 

3c. The Selection Rules 

We shall next consider the laws which govern the interchange of energy 
between an electron and the phonons. There are two important laws, closely 
analogous to the laws of conservation of energy and momentum for two 
masses in collision. The quantity analogous to momentum for the phonon 
is a vector, called P y , directed along the direction of propagation of the 
phonon, and having a magnitude given by the relationship between mo- 
mentum and wavelength 

P y = h(l/\) = h/\, Py || propagation direction (3.3) 

In a transition in which an electron exchanges energy with the phonons, 
and changes its momentum from Pi to P 2 , so that 8n y = ±1 for one of 
the modes, one selection rule requires that 

P» - Pi + StlyPy = 0. (3.4) 

This is analogous to conservation of momentum; actually it is based on far 
more subtle effects. The conservation of energy requires that 

S 2 + hriyhvy = Si (3.5) 

where 

S 2 = Pl/2m, Si = Pl/2m (3.6) 

are the electron's energies before and after collision. The mass m need not 
be the mass of an electron but may instead be the "effective mass," a mass- 
like quantity of the same order as the electron mass which takes into ac- 
count the influence of the periodic potential of the crystal structure upon 
the electron wave packet. 
The effective mass concept represents a simplification that may not 
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necessarily be correct. In a cubic semiconductor, the electron waves can 
te "refracted" as are the longitudinal and transverse acoustical waves. 
The deviations from Ohm's law of Fig. 2 furnish evidence that the simpli- 
fied assumption of equation (3.6) must in fact be replaced by the more 
general possibility. We shall return briefly to this point in Section 5. 

In addition to the conservation of energy and momentum, there are two 
other approximate selection rules which, while not exact, are so nearly 
fulfilled that no appreciable error is introduced by using them: 
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Fig. 4— Comparison of the scattering by acoustical phonons with the scattering of a 
small mass in elastic collision with a larger mass. 

Only Stiy = ±1 is allowed. 

For the acoustical modes, only the longitudinal modes interact with 
electrons. (This restriction does not apply to optical phonons.) 

Figure 4 shows the allowed transitions for an electron with initial mo- 
mentum Pi in the ^-direction. If the energy of the phonons were zero, the 
allowed transitions would be to points on the sphere (or circle in Fig. 4) 
with P 2 = P\ . Since the energy of a phonon is 



kVy = IIC/X = CPy = C \ P 2 — Pi 



(3.7) 



however, the end points lie on the surfaces shown. 

These surfaces do not differ much from the sphere, as may be seen by 
considering the final energy for an electron that reverses its motion by 
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phonon absorption. For this case 

P y = P 2 + Py (3.8) 

and 

(Pi - Pl)/2m = cP y (3.9) 

so that the change in magnitude of momentum is 

Po - P, = 2mc. (3.10) 

For an electron with energy kT and "thermal velocity" 

v T = (2kT/myi\ (3.11) 

corresponding to 10 7 cm/sec at 300°K, the fractional change in momentum is 

(P 2 _ Fl )/ Pl = 2c/v T - 2 X 5.4 X 10 5 /10 7 = 0.11. (3.12) 

Thus the phonon absorption surface lies only 11% outside the constant 
energy sphere. Figure 4 is drawn for the case of P\ = 6 mc, or v x = v T for 
32°K, for purposes of exaggerating the differences in the surfaces. (A fur- 
ther discussion of Fig. 4 is given in the appendices.) 

For transitions with optical phonons, in the range of interest, hv y is 
nearly independent of P y . Furthermore, the optical phonons have hv y = 
k 520°K so that they are nearly unexcited at room temperature and have 
n y = so that only 8n y = +1 is allowed. For this case transitions can 
occur only if 8i is greater than hv oP and the end surface is a sphere with 

S 2 = Si - hv p • (3-13) 

We shall neglect the role of the optical phonons until after a comparison 
between acoustical phonon processes and experiment has been made. We 
shall then show that they play an essential role in explaining Ryder's data. 

3d. Energy Exchange and The Equivalent Sphere Problem 

We shall here give in brief some results derived in the Appendices which 
permit us to show the equivalence of the problem of acoustical phonon scat- 
tering to a problem in gas discharges. This has two advantages: it enables 
us to take over the solution to the statistical problem from gas discharge 
theory, the second problem mentioned at the end of Section 3a, and to 
concentrate on the problem of the mechanism. In addition the equivalence 
makes it much easier to visualize the mechanism of energy losses. 

According to the theory of phonon scattering, an electron is equally likely 
to be scattered from its initial direction of motion to any other. This implies 
that after an interaction the electron is equally likely to end in any unit 
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area of the surfaces of Fig. 4. The probability of being scattered per unit time 
is simply 

1/ti = Vx/t (3.14) 

where Vi is the speed and I the mean free path ; according to the theory I is 
a function of the temperature T of the phonon system and is independent 
of Vi . The time t\ is the mean free time or average time between collisions. 
Figure 4 shows the average energy after collision. The Figure represents 
a case in which the average energy is somewhat smaller than the initial 
energy. This will be the case for a high energy electron, that is one with an 
energy greater than kT for the acoustical modes. The average loss in energy 
for a high energy electron is found to be 

<58> = -c°-P\/2kT (3.15) 

where P y is the momentum change in the collision. This formula is analogous 
to the formula for energy loss if a light mass m strikes a heavy stationary 
mass M and transfers a momentum P*_ — P\ = P y to it. The energy transfer 
is given by (3.15) if 

M = kT/c 2 (3.16) 

since then the kinetic energy of the large mass is simply 

P\[2M = c 2 P\/2kT. (3.17) 

The value of the mass which satisfies equation (3.16) for room temperature 
may be calculated from the previously quoted values of v T and c: 

M = kT/c- = mv%/2c = 170m, (3.18) 

a value which may certainly be considered large compared to m. 

Equation (3.15) is not the complete expression for average energy change 
for a collision with momentum change P y and another term representing 
energy gain also occurs. If the complete expression is averaged over all final 
directions of motion, it is found that the average change of energy, which is 
obviously the average energy change per collision, is 

<58> = 4wc 2 (l - Pl/AmkT) = (4mkT/M) - (Pl/M). (3.19) 

This is the correct expression for the average gain in energy per collision of 
a light mass m colliding with a heavy mass M which is moving with the 
thermal energy appropriate to temperature T. This corresponds to a thermal 
velocity of 

Vtm = (2kT/My- = 2 ,/2 c (3.20) 
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for the large mass. The second term in (3.19) is just the average of (3.15) 
over all directions of motion after collision and represents the energy loss 
that would arise if M were initially stationary. The first term represents an 
energy transfer from M to m due to the thermal motion of the large mass. 

Furthermore, if the light and heavy masses are perfectly elastic spheres, 
the scattering of m will be isotropic, just as is the case for the phonons. This 
shows that there is an almost perfect correspondence between the two 
mechanisms of scattering so that we are justified in using previously derived 
results for the sphere case and applying them to the phonon case. 

To complete the equivalence we should introduce a density of large spheres 
so as to get the correct mean free path. There is nothing unique about this 
procedure, as there is about the mass M and temperature T, and we may 
make a large selection of choices for number of M-spheres per unit volume 
and radii of interaction so as to obtain the desired mean free path. Once 
any choice is made, of course, it will give the same mean free path inde- 
pendent of electron energy and may be held constant independent of the 
electric field. 

3e. Acoustical Phonons and Electric Fields 

We shall next give a very approximate treatment of mobility in low and 
high electric fields. The emphasis will be upon the interplay of the physical 
forces, the mathematical details being left to the Appendices or to references. 

In the Ohm's law range, the field E is so small that the electrons have 
the temperature of the lattice. They have a velocity of motion of approxi- 
mately 

v T = (2kT/m) l i* (3.21) 

and a mean free time between collisions of 

r = 1/vt . (3.22) 

The electric field accelerates the electron at a rate 

a = qE/m (3.23) 

and imparts a velocity ar in one mean free time. Since the collisions are 
spherical, the effect of the field is wiped out after each collision. The drift 
velocity is thus approximately 

Vd = ar= (qC/mv T )E. (3.24) 

An exact treatment which averages over the Maxwellian velocity distribu- 
tion gives a value smaller by 25% and leads to 

Mo - Iqt/Whnvr ■ (3-25) 
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Since theory shows that I varies as T~ x , the mobility should vary as T~ z 2 . 
This prediction is in good agreement with experimental findings over the 
range of conditions for which the dominant scattering processes are those 
considered here. 

Next we consider the effect of very large fields. Under these conditions 
an electron drifting in the direction of the field with drift velocity Va acquires 
energy from the field at an average rate 

(d€,/dlU l0 to a - ViqE. (3.26) 

If this power is large enough, the electrons will be unable to dissipate energy 
sufficiently rapidly to the phonons that they can maintain their normal tem- 
perature. As a result their average energy mounts, after the field is initially 
applied, until they can furnish energy to the phonons fast enough to main- 
tain a steady state. Under these conditions the sum of the two rates is zero 

(</£/<//)duo to B + (.dZ/dl)a U0 to phonona = 0. (3.27) 

If the field is high enough, there may be no steady state solution. This 
can occur if the ability of the phonons to remove energy decreases with 
increasing energy. Such cases play an essential role in the theory of dielectric 
breakdown. 8 In them it is concluded that electrons will gain sufficient energy 
from the field so that they can produce secondary electrons which repeat the 
process thus producing avalanches. For the cases with which we are con- 
cerned, theory indicates that the energy losses increase rapidly with the 
energy of the electron while the power input decreases because of decreasing 
mobility so that a steady state will thus occur. 

In order to estimate the drift velocity for the steady state we must intro- 
duce expressions for the two powers involved. For this purpose we assume 
that an electron has on the average a speed V\ and we calculate the power 
to phonons as the average energy loss per collision for this velocity times 
the rate of collision, V\/ (. For V\ » v T , we can neglect the effect of motion of 
the M spheres and thus obtain from (3.19) 

(rfS/rf/)phonoas = -(»i/0 mhl/M. (3.28) 

The mobility will be less because of the higher collision rate so that the drift 
velocity in the field will be approximately 

v d = (qt/mvi)E. (3.29) 

The power furnished by E will be 

(d&/dl)auo to b - tfe/mvJE*. (3.30) 

8 See the references to Frohlich and Seitz in Section 1 . 
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The steady state condition then leads to 

2-! = {qCE/m) ll \M/m) m (3.31) 

and to 

v d = (qCE/m) l, \m/M) m 

= (VicqCE/mvr) 112 ££ (c/xo£) 1/2 . 

The treatment 9 based on accurate statistics for the equivalent sphere model 
leads to 

v d = 1.23(c M o£) 1/2 . (3.33) 

The transition between the high field behavior and low field behavior 
should occur in the neighborhood of a critical field E c at which both limiting 
forms give the same Vd : 

Vd = noE c = 1.23(cMo£ e ) 1/2 , (3.34) 

leading to 

E c = 1.51 c/mo (3.35) 

and to a drift velocity, which shall be referred to as the critical velocity, of 

v de = 1.51c (3.36) 

if Ohm's law held to a field as high as E c . The drift velocity can be ex- 
pressed in terms of E c by the equation 

v,, = m {EE e )™ (3.37) 

for values of E much greater than E c . 

It is interesting to note that this initial field is just that which would give 
electrons a drift velocity corresponding to the thermal motion of M-masses. 
This seems a natural critical field. For it the effect of random motion of M 
would be suppressed by the systematic drift velocity so that the transfer 
of thermal energy to the electrons would be much reduced. This value of 
E c corresponds to much smaller initial fields than are sometimes proposed. 
For example, one frequently encounters proposals that Ohm's law should 
hold up to the condition that v d = v T . This would correspond to 10 times 
higher field at 300°K than that obtained. Another criterion is that the 
energy gained in one mean free path, qCE, should be equal to kT. This is 
substantially equivalent and corresponds to v d = v T /2. 

B Druyveslcyn Physica 10, 61, 1930. This paper is reviewed by S. Chapman and T. G. 
Cowling in "The Mathematical Theory of Non-Uniform Gases," Cambridge at the Uni- 
versity Press, 1939, page 347. (The factor is 0.897 (18tt/8) 1/4 = 1.23.) 
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For comparison with experiment we note that for a value of 

E = 4E C = 6.04c/mo (3.37) 

such that if Ohm's law held 

v d = 6.04c, (3.38) 

the value of va should be less than half the value predicted by Ohm's law. 
We shall shortly discuss the discrepancy between this prediction and ex- 
periment. 

The "temperature" of the electrons may be conveniently expressed in 
terms of the ratio E/E c . Since the electrons for the high field case are not 
in a Maxwellian distribution of velocities, one cannot define their "tem- 
perature" unambiguously. As a measure of their temperature we shall take 
their average kinetic energy divided by k. This leads to a ratio of electron 
temperature T(E) to crystal temperature T of 2vl/3v T ■ Since to a first 
approximation the ratio of mobilities at low and high fields is 4vi/3ir V2 v T , 
the ratio of temperatures is 

T(E) = 3*[ vol* = 3jE 
T 8 lv d (E)/E] SE e • [6AU) 

This ratio may also be thought of in terms of the square of the ratio of 
drift velocity on the extrapolated Ohm's law line to the drift velocity on the 
E 112 line: 

T(E)/T = (3WmnoE/v d {E)Y. (3.41) 

Either of these equations may be used to estimate electron temperature 
from the data in the range in which the E in formula is a good approximation. 

4. Comparison Between Theory and Experiment 

4a. Discrepancy in Critical Field 

In Fig. 5 we repeat Ryder's data of Fig. 2 together with data on an addi- 
tional sample at 77°K. This new sample is considered more reliable than 
the first since its low field resistivity varies in just the proper ratio [see 
(3.25) and subsequent text] of (298/7 7 ) 3/ ' 2 compared to its valus at room 
temperature. Also we show the theoretical curves that will be discussed 
below. 

The deviations of the data from Ohm's law do not occur at fields as low 
as those predicted in Section 3e. For c — 5.4 X 10 5 cm/sec, the critical 
drift velocity should be 

v dc = 1.51c = 8.2 X 10 5 cm/sec. (4.1) 
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It is seen that Ohm's law is followed to several times higher velocities with 
negligible deviations. The deviations should be a factor of 2 at the field 
corresponding to 



v d = 6.04c = 3.26 X 10 6 cm/sec. 



(4.2) 



on the Ohm's law line. The deviations are actually much less. 

Another important difference between the data and the theory of Section 
3 is that the experimental points do not continue on a straight line with 
slope 1/2 but instead tend to flatten out with a roughly constant drift 
velocity. 
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Fig. 5 — Comparison of E. J. Ryder's experimental data and the statistical theory of 
Appendix A7. 



Three theoretical curves are shown. These are based on an approximate 
treatment that includes the effect of the optical transitions. Due to the ap- 
proximation, the optical modes are neglected below the points marked 
Op on the Figure. This approximation also leads to a discontinuity in slope 
at these points; in a more accurate treatment, this bump would be smoothed 
out. The optical modes play the least role for the curve at 77°K and for this 
theory fits experiment within experimental accuracy if a value of 



Vd, 



= 2.6 X 10 6 cm/sec. 



(4.3) 
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is used. This value is 3.2 times larger than the value given by equation 
(3.36) using c = 5.4 X 10 5 cm/sec, the value appropriate for longitudinal 
phonons. 10 The interpretation of this discrepancy, which we refer to as the 
"low field" discrepancy, is discussed in Section 5. It does not, of course, 
imply an error in the value of the sound velocity, but instead an error in the 
theory leading to the formula for critical velocity in terms of sound velocity. 

Although an exact theory along the lines discussed in Section 5 has not 
been developed, it appears evident that its chief effect will be to increase 
energy interchange with the phonons by a factor of 3.2 squared or approxi- 
mately 10. This increase can be effected in a mathematically equivalent way 
by introducing an effect he Telocity for dealing with phonon energies which is 
3.2 times larger than the true velocity of longitudinal waves. The approxi- 
mate theory in the Appendices uses this procedure. 

We may remark in passing that only two constants were arbitrarily chosen 
to fit the curves to the data. One of these was the effective velocity c = 1.73 
X 10 6 cm/sec. which is 3.2 times larger than the speed of longitudinal waves. 
The other was the mobility of electrons at room temperature. Three other 
constants were chosen from independent estimates of the properties of the 
crystal. One of these is hv for the optical modes for which a value of &520°K 
was used; another is the effective electron mass, for which the free electron 
mass was used; and a third was the interaction constant for optical modes, 
which was set equal to that for the acoustical modes. The meaning of these 
terms is discussed in the Appendices. 

4b. The Ejjccl of the Optical Modes 

We shall next discuss briefly the role of the optical modes before remark- 
ing on a theory of the low field discrepancy. 

As discussed above the optical modes can act only if £ t > hv op . Theory 
indicates, however, that when they do come into action they are much more 
effective than the acoustical modes. On the basis of these ideas, we can see 
how they can act to give a limiting drift velocity that does not increase 
with increasing electric field. For purposes of this illustration we shall imagine 
that so high an electric field is applied that an electron may be accelerated 
from Pi = to P 2 = (2mhi> op ) u -, at which its energy equals hp op , in so short 
a time that it is not scattered by acoustical modes. As soon as it reaches 
Pi , we assume that it is scattered by the optical modes, loses all its energy 
and returns to zero energy. This process then repeats, the period being 

10 This is the velocity of longitudinal vaves in the Til 0] direction as reported by W. L. 
Bond, W. P. Mason, H. J. McSkimin, K. M. 01?en and G. K. Teal, Pliys. Kev. 73. 5W (1948V 
See "Electrons and Holes in Semiconductors," page 528, for the reason for using this wave. 
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P 2 /qE since dPi/dt = qE. The average momentum is evidently P 2 /2 and 
the average velocity is 

v d = P 2 /2m = (W2*»)" 2 (4-4) 

and is independent of E. The optical modes correspond to a temperature of 
about 520°K and this leads to 

v d = 6.3 X 10 6 cm/sec. (4.5) 

in general agreement with the observed value. 

The theory in the appendices indicates that both optical and acoustical 
modes are active simultaneously and their interplay leads to the theoretical 
curves shown. (In the Appendices a further discussion is presented and some 
additional data are compared with theory.) 

The tendency of the theoretical curves to fall below the data for 193°K 
and 298°K for field values below the optical point is thought to arise largely 
from the approximations employed in the theory. The approximations 
neglect the ability of the optical modes to enable the electrons to lose energy 
for fields below the indicated value. Actually some electrons will be scattered 
by the optical modes and this will contribute in an important way to hold- 
ing the temperature down and the mobility up. A correct treatment would, 
therefore, raise the theoretical curve appreciably in the region where it 
deviates most from the data. 

4c. Electron "Temperatures" 

From the theory it follows that the average electron energies correspond 
to about 520°K at the points marked Op on Fig. 5. The highest point on the 
298°K curve corresponds to ~700°K and the highest point on the 77°K 
curve corresponds to 550°K. For this last case the electron temperature is 
more than seven times as high as that of the atomic vibrations. In the ap- 
pendix we quote some other earlier data of Ryder's that indicates electron 
temperatures of about 4000°K while the crystal itself remains at room 
temperature. 

5. An Explanation of the Low Field Discrepancy 

The failure to deviate from Ohm's law at the low fields predicted indicates 
that the electrons can dissipate their excess energy more effectively than 
would be expected on the basis of their mobility. This conclusion is forced 
on us by the observation that they apparently retain their thermal dis- 
tribution and normal mobility to higher fields than predicted. It is not pos- 
sible to explain the discrepancy by assuming a large or a small value for the 
effective mass, since the value of the effective mass does not enter into the 
final comparison with experiment. 
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It is possible, however, to explain the discrepancy by assuming that the 
effective mass is not single valued. This assumption corresponds to the case 
in which the surface in the Brillouin zone belonging to a single energy is not 
a sphere but instead a complex surface of two or three sheets. Such surfaces 
have been found as a result of numerical calculations for certain crystals 11 
and it has also been shown that such surfaces are to be expected in general 12 
if the energy at the bottom of the conduction band is degenerate. It appears 
necessary to assume that such complex surfaces occur in order to explain 
magnetoresistance effects. 13 

In terms of Fig. 4, this theory replaces the circular energy contours by 
deeply re-entrant curves. Transitions from peak to peak of the curves result 
in large energy transfers to the phonons and hence more effective energy 
losses. This effect can occur without a compensating change in the effective 
mass involved in the mobility and, as a result, the critical field may be 
increased by a large factor. A preliminary analysis indicates that in order 
to increase the critical field by a factor of 3 a value of about 3 is also 
required for the ratio of maximum to minimum momentum for the energy 
sy^face. A similar analysis of magnetoresistance leads to a factor about 50% 
larger in order to account for the increases in transverse resistance of about 
7-fold observed by Suhl. 14 At the time of writing, therefore, the author feels 
that both the critical field data at low fields and the magnetoresistance data 
require a modification of the effective mass picture and that the same modi- 
fication may well explain both sets of data. 

I am indebted to E. J. Ryder, whose experimental results provoked the 
analysis presented in this paper, to F. Seitz and J. Bardeen for several helpful 
discussions, to Gregory Wannier for an introduction to the analogous case 
in gas discharge theory and to Esther Conwell for help with the manuscript. 

I shall also take this opportunity to express my appreciation to C. J. 
Davisson. The opportunity to work in his group was a large factor in my 
decision to come to Bell Telephone Laboratories, where I enjoyed his stimu- 
lating companionship while assigned to his group, and later as well. 

APPENDICES 

A.l Introduction and Notation 

The problem of energy exchange between the electrons and the phonons 
recjuires a somewhat more sophisticated treatment than does the problem 
of mobility at low fields. In order to present the theory of energy exchange, 

II W. Shocklev, Pliys. Rev. 50, 754 (1936). 

12 W. Shockley, Pliys. Rev. 78, 173 (1950). 

13 W. Shockley, Pliys. Rev. 79, 191 (1950). 

14 H. Suhl, Pliys. Rev. 78, 646 (1950). Suhl finds increases in resistance in transverse 
fields as high as 7-fold. 
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it is necessary to reproduce a large amount of the material dealt with in or- 
dinary conductivity theory. We do this in a somewhat abbreviated form ex- 
panding the exposition on the points particularly pertinent to the theory 
of energy losses. 

For convenience we reproduce here a number of the more important sym- 
bols. The references indicated refer to places where they are discussed in the 
text. 

a = lattice constant; Fig. 3. 

c = speed of longitudinal acoustical wave; Equ. (3.2). 
en = average longitudinal elastic constant; Equ. (A4.1). 

e = base of Naperian logarithms. 
E = electric field. 

S = energy. 

Si„ and £ 2n ; Equ. (A4.1) and (A7.9). 

h = 2tt h = Planck's constant. 

k = Boltzmann's constant. 

I = mean free path for electron due to scattering by acoustic phonons; 
(A4.3). 
/ op = describes scattering by optical phonons; (A7.19). 
m = effective mass of electron. 
M = mass in equivalent mass treatment; (A5.8). 

P = "crystal momentum" of electron = // times its wave number. 

V = volume of crystal. 

A = dilation; (A4.1) and (A7.10). 

v = frequency of normal mode. 
v op = frequency of optical mode (used in Section 4 only); Equ. (4.5). 

A .2 The Probability or Transition into Energy Range 58 2 

In this section we consider an electron initially with energy S x and mo- 
mentum P u which for convenience we take to be along the iVaxis, and 
we evaluate the probability that it make a transition to states with ener- 
gies in the range S 2 to 8 2 + 5S 2 . We shall assume that the crystal is elas- 
tically isotropic so that for the spherical energy surface approximation 
employed, i.e. equation 3.6, the scattering will be symmetrical about the 
iVaxis. The end states, K, may, therefore, be consideredm groups lying 
in the range </S 2 , dd where 6 is the angle between P x and P v These states 
lie in a ring in JF*-space whose volume is 

2irP 2 sin P 2 dd dP 2 = 2wtnP 2 sin d dd JS 2 (A2.1) 
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The number of end states in dd d& 2 space is thus 15 

(V/h 3 ) ImnPt sin 6 dd dZ 2 = P dd dZ 2 (A2.2) 

(The density p introduced above is used below in calculating the transition 
probability; since spin is conserved in the transitions of interest, the den- 
sity of possible end states in phase space is \/h instead of 2/h .) 

The transitions will occur between states of the entire system, electron 
plus phonons, which conserve energy. The transition of the electron from 
Pi to P 2 requires a compensating change in the phonon field. 16 The con- 
servation laws allow two possibilities: (I) phonon emission; the longitudinal 
acoustical mode with 

P a = hka = -(Pa- Pi) (A2.3) 

undergoes a change 

«« -> »« + 1 (A2.4) 

with a change in energy for the electron of 

82 - Si = -h<a a = -hc/\ = -cP a (A2.5) 

where c is the velocity of the longitudinal phonons that are chiefly respon- 
sible for the scattering. These relationships lead to conservation of the 
sum of P for the electron plus /2 n a P a for the phonons. The other pos- 
sibility is (II) phonon absorption, for this case 

f e m hkf, = (P 2 - K) (A2.6) 

iip ->n - 1 (A2.7) 

82 - 81 == +cP (A2.8) 

again with conservation of the sum of P vectors. 

If we denote by £ the energy of the electron after collision plus the change 
in phonon energy, then the requirement of equality for energy before and 
after collision gives 

& = & 2 + SllyCPy = Si (A2.9) 

where 8n y = +1 is the phonon emission or a surface and 8n y = — 1 is the 
phonon absorption or /3 surface of Fig. 4. 

15 The notation in this appendix follows closely that of \V. Shockley, "Electrons and 
Holes in Semiconductors," D. van Noslrand (1950) to which we shall refer as E and II 
in S. See page 253 for a similar treatment of p. 

10 This condition is analogous to one for the conservation of momentum but has a 
different interpretation. See for example E and II in S, p. 519 equation (15). 
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We shall next insert these symbols into the conventional expression for 
transition probability. We consider a system described by one or more 
sets of quantum numbers, say Xi , x 2 , ■ ■ ■ , x n which may take on discrete 
but closely spaced values so that the numbers of states lying in a range 
dx\ , ■ ■ ■ , dx n is 

p(x u ■■■ , x n ) d Xl ■ ■ ■ dx n - (A2.10) 

The system may make a transition from an initial state <po and energy So 
to another state ipi of the same energy between which there is a matrix 
element £/„,-. The total probability of the system making a transition per 
unit time to the range of quantum numbers dx 2 , dx 3 , ■ ■ • , dx„ is then 

W oi dxt ■ ■ ■ dx„ = (2ir/h) | U 0i \- [p/id&i/dxi)} dx 2 ,■■■ , dx n (A2.ll) 

where d&i/dxi is evaluated where S,- = 8 ; if for the range dx 2 , ■ ■ ■ , dx n of 
the other quantum numbers there is no X\ value that gives £,- = So, then 
the transition does not occur. 

We shall apply this to our case letting 9 = *i and S 2 = *2- The expres- 
sion dS,/d.Vi then becomes 

*Lm(^\ =8n y c( d -^) (A2.12) 

ax t \ddj Si v do y g2 

where 

£ * - £ l R - * l - 1 crt + p\ - 2p.p, cos ey- (A2 13) 

= P } Pi sin d/Py . 

We then obtain, for W n d& 2 , the probability per unit time of transition of 
the electron from I\ to states with energies between S 2 and S 2 + <&2 , the 
expression 

W 12 dS, 2 = (2ir/h) | U | 2 (V/h 3 )2irmP 2 sm 6 X (P y /c8n y PiP2 sin 6) d& 2 

= (V/2irh 4 )m | U | 2 (P y /cSny Pi) dS 2 (A2.14) 

= (V/2-nh 4 )m | U | 2 (Py/PJi-dPy); 

where the negative coefficient of dP y is without significance except for its 
relationship to the selection of the limits of integration. In subsequent equa- 
tions we shall disregard the sign convention which relates d& 2 to dP y ; no 

"See L. I. Schiff "Quantum Mechanics," McGraw-Hill Book Co. (1949), equation 
(29.12). The additional factor L/(3£«'/<toi) converts the p used here to that of Schiff, winch 
latter is number of states per unit energy range. 



HOT ELECTRONS IN GERMANIUM AND OHM'S LAW 1015 

error is introduced provided the subsequent integrations are always in the 
direction of increasing values for the variables concerned. 

A. 3 The Allowed Ranges for P y 

An electron with an energy corresponding to room temperature can 
change its energy by only a small fraction in a one phonon transition. The 
extremes occur for d = ir corresponding to complete reversal of direction. 
For this case we have 

S 2 - 8, = (Pa - P\)/2m = -8)iycP y 

(A3.1) 
= -bn y c{P«+ PO 

so that 

Pi - Py = - bn y 2mc = - 25» T P . (A3.2) 

Thus the limiting values of P 2 differ from Pi by 

±2P = ± 2mc (A3 .3) 

in keeping with the results shown in Fig. 4. [For Vi = P\/m = 10 cm/ sec, 
corresponding to 81 = 0.025 electron volts, and c = 5.4 X 10 5 cm/sec, it is 
seen that P 2 and P\ differ by 10%.] For this case the range of P y is 

phonon emission, 8>i„ = +1, P. from to 2 (Pi - P ) (A3.4) 

phonon absorption, 8ng = — 1, Pg from to 2(Pi + P a ). (A3. 5) 

A singularity occurs for Pi = P . For this case phonon emission becomes 
impossible and the inner curve of Fig. 4 shrinks to zero; in Fig. Al we show 
the sequence of shrinkage. The value of 81 for this condition corresponds to 
thermal energy for a temperature of less than 1°K. Under the conditions 
for which we shall compare theory and experiment, a negligible number of 
electrons lie in this range. Accordingly we shall use the above limits in cal- 
culations and neglect the small errors introduced. 

A. 4 The Matrix Element and the Mean Free Path 

The matrix element may be written in the form 

I U | 2 = 8L (2n y + 1 + 8ny)Pyc/4Vcu (A4.1) 

where 8i„ is the derivative of the edge of the conduction band in respect 
to dilatation of the crystal and Cu is the elastic constant for longitudinal 

18 See E and H in S, page 528, equation 31. The second expression in equation 31 
of this reference is in error by omission of a factor hoika = cP y . 
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Fig. Al -Initial and Final Values of P for Transitions which conserve energy. 

(a) The two nearly spherical surfaces for | Pi | =6 mc. 

(b) The two surfaces for |jP, | = 3 mc. 

(c) For an electron with | K I = mc, energy loss is impossible. 

(d) The case of | F, | = 0. 

waves. Inserting this in the expression (A2.14) for Wu dS 2 , we obtain 
W12 dZi = (F/2ttM) m [e,l n P y c/iVc(t] 

X (2n y + 1 + 8ii y )(P y /Pi) dP y (A4.2) 
= (l/0(l/8mP 1 )(2« 7 + 1 + 8, h )(cP y /kT)P y dP y 
where we have used the symbol I to represent 

I = T h A c*//m 9 8l, kT (A4.3) 

because, as we shall shortly show, I is the mean free path for electrons. 
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For the cases with which we shall be concerned, the values of n y may 
be approximated by classical equipartition. This may be seen from the 
fact that largest energy phonons correspond approximately to an energy of 

cP y = 2cPi = (icm/Pi)Pl/2m 

(A4.4) 
= (4e/vi)kT. 

Their energies will, therefore, be considerably less than kT. For large in- 
creases in electron energy in high fields, however, this approximation may 
not be adequate. At room temperature cP y /kT = 0.2 and the critical 
range will correspond approximately to an increase of about 10 fold of 
electron temperature above the temperature of the crystal. Under these 
conditions cP y deduced from equation (A4.4) will be about 2kT for the 
most energetic phonons; for this condition, however, P y lies at the edge of 
the Brillouin zone and dispersive effects must be considered. In this treat- 
ment we shall not investigate further these limits and shall in general 
assume that cP y < kT. 

We shall next derive an expression for the mean free path and verify 
that the scattering is isotropic. These results can be derived more simply 
and directly from the matrix element by neglecting (P /Pi) and (cP a /kT) 
from the outset. For the treatment of energy losses that follows, we cannot 
make these approximations. We shall, however, make them in the re- 
mainder of this section thus establishing that our more general formulation 
reduces correctly to the more convenient and simpler formulation usually 
used. 

For the condition under which equilibrium applies we may approximate 
n y as follows: 

n y = l/[(exp cP y /kT) -1) = (kT/cP y ) - \. (A4.5) 

Then, for the phonon emission or a case, the contribution to Wu d&z be- 
comes 

W n d& 2 = {±mtPi)- l {\ + (cP a /2kT)]P a dP a (A4.6a) 

and for the phonon absorption case, it becomes 

WndZ, = (4w/'P 1 )- , [l - (cP l3 /2kT)]Pf 3 dP t) . (A4.6b) 

We shall use these expressions later for the calculation of energy exchange, 
in which case the terms in cP a /2kT, which favor phonon emission, play an 
important role. In order to check the expression for mean free path we neg- 
lect these terms, however, and also approximate the integral of P a dP a by 
2P\ rather than 2 (Pi — P ) 2 . The total probability of transition from state 
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Pi, which should be taken to be l/n where n is the mean free time, is 
then 

l/ ri = Wi= (AmlPv)- 1 4Fi = (Pi/m)/( = vJL (A4.7) 

This is just the relationship appropriate to the interpretation of I as a mean 
free path between collisions. 19 It does not follow that ti is the relaxation 
time for the current, however, unless the average velocity after collision 
is zero. We shall next show that the average velocity after collision is zero 
to the same degree of approximation used above by showing that scatter- 
ing into any solid angle of directions is simply proportional to the solid 
angle, i.e. the direction of motion after collision is random. 

The conclusion that the scattering is nearly isotropic follows from the 
approximate proportionality of probability to P y dP y . Since P 2 is nearly equal 
to Pi and substantially independent of 0, we may write 

P a dP a = \ d(P a Y = h d (2P\ - 2P\ cos 0) 

= - P\d cos 6 = P\ sin dO. (A4.8) 

The last term is simply proportional to the solid angle lying in range d6; 
hence the end states are distributed with uniform probability over all direc- 
tions and the scattering is isotropic. 

A5. Approximate Equivalence to Elastic Sphere Model 
In this section we shall show that on the average the energy exchange 
between the electron and the phonons when the electron is scattered through 
an angle is very similar in form to that corresponding to elastic collisions 
between spheres with the phonons represented by a mass much greater 
than the electrons. If dP a and dPp correspond to scattering through angles 
between and + dd, then the energy loss for phonon emission is cP a and 
the energy gain for absorption is cP p . The relative probabilities of loss and 
gain are given by equations (A4.6) and from these it is found that the 
average energy gain is 

, v cP p [l - (cP /2kT)]P g dP - cP a [l + (cP a /2kT)]P a dP a 
{88) " [1 - (cP /2kT))P dP, + [1 + (cP a /2kT)]P a dP a 

_ c[Pl dP,, - P\ dP a \ - (c 2 /2fer)(PJ dP + Pi dP a ) 
~ [1 - (cP /2kT)]P dP B + [1 + (cP a /2kT)]P a dP a ' 

Since we are concerned chiefly with cases in which Pa « Pi, Pi — P\ , 
and cPft/kT « 1, we may set 

P a = P S = /\ = 2P, sin (0/2) (A5.2) 

19 See E and H in S, Chapter 11. 
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where P y corresponds to neglecting the energy change of the electron on 
collision. This approximation is not good enough for the first term in (A5.1) 
which involves Pp — P a . In order to evaluate this first term we note that 
the relationships 

P a = P y [l - (P P y /2Pl)] (A5.3a) 

Pfl = Py[l + (P Py/2P\)} (A5.3b) 

may be derived up to the first order in P . This permits us to write 

,2 ,r, n 2 ,n . /i x , r/ r. i. W-> n2 



(A5.4) 



pi dp e - pi dp a ± (|) d [(p, - p„)(3p* t >] 

= rf[PoFi/PS] - 4(P P 2 7 /-Pi) /\ <^V 

Hence 

<5S) = 2(cP P 2 7 /-Pi) - c 2 P\/2kT. (A5.5) 

The second term is proportional to the (average change in momentum) 2 
for collision by angle 6. It thus corresponds to energy which would be trans- 
ferred to an initially stationary mass M » m by the colliding electron pro- 
vided we take 

M = kT/c\ (A5.6) 

That this mass is much greater than the electron's mass may be seen in 
terms of v, the velocity of a thermal electron: 

mv\/2 = kT. (A5.7) 

From this we obtain 

M = w(» 1 /c) 2 /2 = 170w (A5.8) 

at room temperature where V\ — 10 cm/sec while c = 5.4 X 10 cm/sec. 
The first term may then be interpreted as follows: 

2cP P\/P\ = 2c"-mP- y /P\ = 2(kTm/M)(P\/P\) (A5.9) 

If this is averaged over all angles 6, the Py/Pi term becomes 2; the energy 
gain is then just that picked up by a mass m colliding with a mass M moving 
with a Maxwellian distribution at temperature T as may be seen as fol- 
lows: For this case the velocity Vu of M parallel to the line of centers on 
collision imparts an added velocity 2v M to the electron and, on the average, 
an energy 

(£)ro <(2z) A ,) 2 ) = 2m (v 2 M ) = 2 mkT/M (A5.ll) 
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since M(i^,)/2 = kT/2. In addition, however, there is an effect due to 
relative velocity: when v M is directed so as to increase the closing velocity, 
the probability of collision is increased. Due to this effect, collisions with 
higher relative velocity are favored and as a result the energy transferred 
due to the v u effect is just doubled'" leading to a total contribution of 

(energy transfer due to v M ) = 4 tnkT/M. (A5.10) 

This is just to the first term of (A5.5) when averaged over all values of 0. 

Thus the average of the gain in energy term in (A5.5) is just that corre- 
sponding to interactions with heavy masses M in thermal agitation. The 
difference is that in the sphere model the average energy gain term is inde- 
pendent of 0, whereas in the phonon case it varies as P\/P\ and approaches 
zero for forward scattering so that the dependence upon angle is different. 
The energy loss term, however, is correctly represented by the sphere model. 

The average value of (5S) averaged over all values of is denoted by 
(5S)/>,. Since P\ averaged over is 2P\, we obtain 

(8£,) Pl = 4cP - c 2 P\/kT 

(A^.12) 

- Ante 2 (1 - Pj/4 mkT). 

From this expression it is seen that an electron with energy P\/2m = 2kT 
keeps the same energy on the average after M collision. We shall use ex- 
pression (A5.12) in Section A.6. 

For high electric fields, the electron energies are higher than thermal and 
the loss terms predominate. Furthermore, the scattering in both cases is 
nearly isotropic if M » m and the colliding particles are spheres. Hence, 
the analysis of kinetic theory of ionized gases can be applied to a high 
degree of approximation to estimate electron behaviors. 

It should be stressed that several approximations are involved in this 
treatment. In particular it is assumed that (I) cP y < kT and that (II) 
Pi » P . If this is true, then 

Vi/c = Pi/Po»l (AS.13) 

so that the approximation used in considering the heavy spheres to be 
moving slowly holds and the mass of the heavy spheres is much greater 
than the electron mass: 

M/tn = (kT/c 2 )/(2kT/v\) = v\/2c. (A5.14) 

If conditions (I) and (II) are not satisfied, the approximations leading to 
(A5.5) will require revision. 

20 If v in is the velocity towards centers, then the probability of collision is weighted by 
[1 + (vu/vi n )] and the term linear in v u in the energy, which is (i»i)(4vu«if)> con- 
tributes 2(v 2 u )m to the average transfer. 
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A. 6 Approximate Treatments of Mobility in High Fields 

A correct treatment of mobility in high electric fields E is based upon 
finding the steady state distribution function /(P, E, T) which satisfies the 
Boltzmann equation, i.e. a function for which the rate of change due to 
acceleration by E just balances that due to scattering. This method leads 
at once to rather formidable mathematics which may tend to obscure 
somewhat the physical forces at work. We shall in this section derive re- 
lationships between drift velocity and E on the basis of simpler models 
and shall compare the results with the exact treatment as given in the 
literature for the case of a gas. 

For this purpose, we shall first suppose that the field in effect raises the 
electrons to a temperature T e which is greater than the temperature T of the 
phonon distribution. The electrons with temperature T e will have colli- 
sions at a rate {T e /T) xl - greater than before and their mobility will be re- 
duced from its equilibrium value no to a new value n 

M = (r/r.) l ' 8 Mo. (A6.i) 

The average rate at which the electric field does work on an electron is then 

(dZ/dt) IMd = Force X Speed = qi iEr. (A6.2) 

For steady state conditions this must be equal to minus the average rate 
at which an electron gains energy from the phonons. Denoting this by 
(</£/<//) phono™ we have 

(</S/rf/)fi„i d + (rfS/<//)„ honon8 = 0. (A6.3) 

In order to calculate the average rate of energy loss to the phonons, we 
consider the average energy gain of an electron of momentum Pi. As given 
by(A5.12): 

(8Z) Pl = 4wc 2 [1 - (vtvl/UT)]. (A6.4) 

According to our assumption, the number of electrons in the velocity range 
v to v + dv is N(v) dv = A exp (—mv 2 /2kT e )v-dv. These electrons suffer 
collision at a rate v/l. Hence the average rate of energy gain is 



(A6.5) 



(</S/<Z/) phonon8 = f (8&) P (v/t)N{v) dv/ Ji\(v) dv 

= (8/vW«*VWl - (v e /v T ) 2 } 
where we have introduced 

v T = (2kT/m) 1 ' 2 , v. = (2kT e /m) 1 ' 2 . (A6.6) 
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We see that for thermal equilibrium, with T = T c and v = v e , this equation . 
gives correctly the result that there is no net interchange of energy be- 
tween electrons and phonons. 
The equation for mobility for the case of phonon scattering is 

Mo = 4^/3 "\/W». (A6.7) 

This expression may be used to reduce the steady state equation: 
= (de,/dt) ne id — (iS/rfOphoMons . 

= {v T /v e )qn*& + (S/Vr)(mc 2 v T /C)(v c /v T ) [1 - {v e /v T y\ 
to 

{vjvrf - {v e /v T f = (3tt/32)(h E/cJ\ (A6.9) 

This equation may be solved for v e /v T : 

{vJv T = V(l)(l + [1 + (3t/8)(moEA) 2 ] 1/2 ) 1 ' 2 (A6.10) 

The drift velocity then becomes 

v d = y.E = M o£ V2/(l + [1 + (37t/8)(mo£A) 2 ]" 2 ) i/2 . (A6.ll) 

For E « c/iio, we have 

v d = n,E. (A6.12) 

For E » c/no, we have 21 

v d = (32/37T) 1 / 4 (mo-Ec) 1 ' 2 = 1.36( M o£c)" 2 . (A6.13) 

These equations define a critical field E at which the two limiting cases 
would give the same mobility: 

E c = (32/3tt)"Vmo = 1.84c/ M o. (A6.14) 

At this critical field, the drift velocity is less than the value on either limit- 
ing form by a factor of 

v d /»oE c = V2 / [1 + 5 1 ' 2 ] 1 ' 2 = 0.785. (A6.15) 

This reduction in mobility corresponds to an electron temperature of 

T e = r/(0.785) 2 = 1.627\ (A6.16) 

Expressed in terms of c the drift velocity given by the limiting forms is 

v c = v, t (extrapolated) = 1.84c (A6.17) 

21 This relationship has been published in Bulletin of Am. Phys. Soc, Vol. 26, No. 1, 
paper 55. 
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and the actual drift velocity is 

v d (E c ) = 0.785 X 1.84c = 1.45c. (A6.18) 

If we introduce E c into the steady state equation (A6.9), and express 
Ve/v T as a mobility ratio, we obtain 

(mo/m) 4 - (mo/m) 2 = (E/E c )\ (A6.19) 

We give this equation in order to show the similarity of the Maxwellian 
distribution case to the cruder case considered next. 

A similar treatment may be given for a hypothetical distribution of elec- 
trons such that all have the same energy 8 and speed v. The effective mo- 
bility of such a distribution is 22 



(A6.20) 



M = (qr/m) (1 + (I) dfnr/dfnv) 

= 2qC/3vm = w 112 /j. v t /2v. 

The steady state equation deduced from (A6.19) and (A6.4) is 

(v/2v%) [(t>V2»r) - 1] = (3tt/64) (ji E/c) 2 (A6.21) 

For high fields we find that this steady state condition gives 

v„ = (tt/3) 1 /' (c M0 £)" 2 = 1.01 (c M0 £) I/2 , (A6.22) 

a value somewhat smaller than that obtained from the Maxwellian approxi- 
mation. This leads to a critical field of 

E c = (tt/3) 1 ' 2 c/mo = 1.03 c/mo (A.623) 

at which v d given by (A6.22) extrapolates to give the same value as n E. 
We may use this distribution and make it give identical results with the 
Maxwellian distribution. If we use the steady state condition for low fields, 
we find v = 2 1/2 v T and 

M = Mo = V^r78 mo = 0.625 mo. (A6.24) 

In terms of this fi , the steady state equation becomes 

(mo/m) 2 [(mo/m) 2 - 1] = (E/E' c y. (A6.25) 

with 

E' c = (8/3) »' 2 c/mo = 1.63 c/mo- (A6.26) 

It is evident that if we chose modified values of c and I so as to make 
Mo(0 become equal to mo and E c (c , t ' ) = E c (c, (), then the monoenergetic 

- See E and II in S problem 8 page 293. 
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approximation will give the same results as the Maxwellian distribution. 
We use this procedure in the following section. 

As pointed out in Section A.5, the problem treated here is closely analo- 
gous to electronic conduction in gasses. For this case, an exact treatment 
has been given for high fields such that the motion of the large masses M 
may be neglected. The drift velocity is found to be" 

v d = {ifi^/^im/Myi^qEC/myiyYH). (A6.27) 

Substituting kT/c 2 for M and using (A6.7) to eliminate t, we obtain 

•■W^* (A6.28) r 

= 1.23 faoEcyi*. 

This value lies intermediate between the two simple approximations con- 
sidered above and leads to a critical field of 

E c = 1.51 c/n (A6.29) 

at a velocity of 

v Je = noE c = 1.51 c (A6.30) 

Since the exact case gives v d = fx E for low fields and v d = no(EE c ) u - for 
large fields, it is evident that either the Maxwellian or single energy dis- 
tribution will approximate it well (provided suitable choices of no or fi and 
E c or E' c are made) except for a small error near E e . 

The distribution in energy for the gas case leads to a probability of 
finding the electron in a range v to v + dv proportional to 

[exp - f tnv dv/(kT + 3M(qCE/3mv) 2 ) v dv (A6.31) 

For high fields this reduces to 

[exp - (Sm/AMXv-m/qCE)-]^ dv 

It is seen that this distribution weights the low energies less heavily than 
does the Maxwellian for the same average energy and thus gives a lower 
mobility. It weights them more heavily than does the single energy and, 
therefore, gives a higher mobility than it. 

23 Druyoestevn Physica 10, 61 (1930). See also S. Chapman and T. G. Cowling, "The 
Mathematical'Theory of Non-Uniform Gases," Cambridge at the University Press, 1939, 
page 351. 
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A. 7 The Effect of the Optical Modes 

A. 7a. Introduction 

The treatment presented above is based entirely upon interaction with 
the longitudinal acoustical modes of the crystal. Since the diamond struc- 
ture has two atoms per unit cell, it is also possible to have "optical modes" 
in which the two atoms vibrate in opposite directions. Such modes may 
have long wave lengths; for example, do the same thing in each unit cell, 
and thus correspond to small values of P y . Hence they may interact with 
the electron waves with P y — Pi. Their frequencies correspond roughly to 
a wave length of about (2) the lattice constant in the flOO] direction and 
hence to a frequency of about 

(4.92 X 10 5 )((i)5.6 X 10- 8 )- 1 (2/t) = 1.12 X 10 13 see" 1 (A7.1) 

The last factor of (2/7r) is a crude allowance for dispersion, which leads to 
a decreasing phase velocity at short wave lengths. This corresponds to an 
energy 

hv = k 520°K. (A7.2) 

These optical phonons thus contain much more energy than the acousti- 
cal phonons; the latter having at room temperature an energy of about 

cP x = 2(cM) (wiPi/2) = (1/10) £ 300° = k 30° K (A7.3) 

[or about k 100°i£ if we use the higher effective value of c discussed in 
Section 4]. Furthermore, the optical phonons will be only slightly excited; 
thus collisions with them will in general involve phonon emission so that a 
collision will on the average result in an energy loss of nearly 500&. This is 
very large compared, for example, to the average loss of about 8 mc 2 — k 
12° per collision for electrons with energies of 4kT. 

A difficulty with the optical modes is that for the approximation of 
spherical energy bands, which we have used in earlier parts of this paper, 
they should have matrix elements which vanish when Pi = 0. This con- 
clusion is reached by considering the deformation potentials corresponding 
to an optical displacement: this may be thought of as moving one of the 
face-centered cubic sublattices of the diamond structure in respect to the 
other. Since the initial position is one of tetrahedral symmetry, there can 
be no first order change in the energy £ c at the bottom of the conduction 
band. There may be a distortion of the energy spheres for higher energies, 
however, and this can lead to matrix elements proportional to P\ and transi- 
tion probabilities proportional to Pi . 

25 This view is in disagreement with the position stated by F. Seitz in his two papers 
on mobility. Phys. Rev. 73, 550 (1948) and 76, 1376 (1949). See for example the text 
between equations (14) and (15) of the latter paper. 
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On the other hand, if the band is degenerate and has energy surfaces 
consisting, for example, of three sheets, then an optical displacement may 
split the degeneracy and the shift in energy for Pi = can be linear in the 
displacement. (For example consider one wave function with angular de- 
pendence of the form (cos 4- sin cos <p + sin sin ip) i.e. a />-type wave 
function with its axis along [111]. Its energy should certainly change for 
relative displacements of the two sublattices along the [111] direction since 
positive and negative displacements are unsymmetrical in their distortion 
in respect to this line.) 

Chiefly from evidence on magneto-resistance, the writer is convinced 
that the electron energy band is complex in form. Thus it would be expected 
that the optical modes would have matrix elements for low values of Pi. 
We shall assume that this is the case but shall not endeavor to deal with a 
non-spherical band. This an inherently inconsistent approach, but should 
give at least a semiquantitative agreement with an exact theory. 

In order to illustrate the effect of the optical modes, we shall assume for 
the moment that for high fields they are the dominant mechanism of scat- 
tering and that the scattering is isotropic. If the mean free time between 
collisions is t, then the power input is 

(d&/dl) fie i d = q*T&/m. (A7.4) 

If the temperature is so low that the optical modes are only slightly ex- 
cited, the transitions will in general absorb an optical phonon so that 

(<fS/rf/) p = hv/r. (A7.5) 

The steady state condition leads to the surprising but simple result that 

v lt = qrE/m = (Iiv/m) 1 ' 2 (A7.6) 

so that the drift velocity is independent of E. 

If we insert k 520°K for hv and the free electron mass for m, v d becomes the 
velocity of an electron with k 260°K of energy giving 

v d = 0.88 X 10 7 cm/sec. (A7.7) 

The limiting value on Fig. 2 for 29S°K corresponds to extrapolating Ohm's 
law to about 1500 volts/cm or a drift velocity of 1500 X 3600 = 0.54 X 10 7 . 
The limiting value for 77°K is about twice as high. These values are seen 
to be in reasonable agreement with the predicted value. 

It should be pointed out, however, that the answer obtained for v d de- 
pends implicitly on the assumption that a relaxation time may be used in 
the simple way employed above. To illustrate that the result is not com- 
pletely general we refer the reader to equation (4.4) which was obtained 
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on the basis of a somewhat different treatment. According to (4.4) 

v d = (i) (2hu/my 2 = (hu/2my' 2 

. (A7.8) 

= 0.63 X 10 cm/sec. 

a result smaller than (A7.6) by a factor of 2 I/2 . 

A correct treatment of the optical modes together with acoustical modes 
would involve solving the Boltzmann equation to find the steady state 
distribution. This will obviously present problems of considerable com- 
plexity. In particular scattering will suddenly begin to increase when the 
electron acquires an energy, 8i > hv and it seems unlikely that analytic 
solutions can be obtained. Even the Maxwellian distribution leads to 
somewhat complicated integrals. 

A. 7b. Estimate of the M atrix Element 

In order to proceed further we must estimate the order of magnitude of 
the optical scattering matrix element. For this purpose we introduce a 
"deformation potential" coefficient for the optical modes by the equation 

fi 2n = dS c /d(x/x ) (A7.9) 

where x is the displacement parallel to the x-axis of one sublattice in respect 
to the other and x Q is the x-component of relative displacement of the sub- 
lattices for equilibrium conditions. The same reasoning as used in treating 
mobility by deformation potentials" may then be applied and the matrix 
element evaluated by analogy with the dilatation waves. For the latter the 
matrix element may be written in the form 

\Ul \- = sL<A 2 )/2 (A7.10) 

where A is the dilatation and (A 2 ) is the average (dilatation) 2 for the 
mode before and after transition." Since half the energy in a running wave 
is potential 

i c„(A 2 ) V = \ hv X [average of (» + J)] 

(A7.ll) 
= hv (2m + 1 + «»)/4. 

This leads to the form introduced in equation (A4.1). By analogy, for the 
optical modes we should take 

I U op | 2 - 8L((V*o) 2 >/2 (A7.12) 

26 VV. Shockley and J. Bardeen Phys. Rev. 77, 407-408 (1950) and J. Bardeen and 
VV. Shockley, Phys. Rev. SO, 72 (1950).' 

27 See E and H in S, page 528. 
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where the stored energy for deformation (x/x ) is 

| coo (x/x y V (A7.13) 

and the average value for a transition from n = to n = 1 is hv for the 
total energy and (\)hv for potential. This leads to 

| t/ p | 2 = fSLkf/lemV (A7.14) 

As a first approximation we may take the stiffness between the planes of 
atoms separated by x as the same as the macroscopic value. This leads to 

coo = c (l (A7.15) 

Furthermore, equal relative displacements of neighbors are produced by 
equal values of A and x/x . Hence approximately equal changes in energy 
may occur so that we may assume that 

S 2n * S ln . (A7.16) 

Under these conditions 

| U op | 2 = (hv/kT) | U A | 2 . (A7.17) 

Since the energy of the optical modes is a maximum for P y = 0, it will 
change only a small fraction for values of P y comparable to Pi. 23 (See Fig. 3.) 
Consequently, conservation of energy leads to transitions between Pi and 
a sphere with P 2 = [2»»(8j - hv)] 112 . The probability is equal to each point 
on the sphere and the transition probability is readily found" to be 

1/rop = (V | U op \*tn*/irh*)Vi (A7.18) 

where v 2 = Pi/m is the speed after collision. If we assume relationship 
(A7.17) between matrix elements and introduce I as defined in (A4.3), then 

l/r p = (hv/kT) v,/Z = v 2 /C op (A7.19) 

where 4 P is a sort of mean free path for optical scattering. 

The dependence of v 2 upon the velocity before collision v x is obtained as 
follows: 

8 2 = Si - hv (A7.20) 

* - [2(51 - hv)/m]m (A7.21) 

= (A - vlyi* 

a See F. Seitz, "Modern Theory of Solids," McGraw-Hill Book Co., 1940, p. 122. 
29 See E and H in S of A.2, p. 493, for a similar treatment. 
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where v, is the velocity corresponding to hv. 

v v = (2hu/my 2 . (A7.22) 

The rate of energy loss to the optical modes is simply 

Wop. (A7.23) 

A. 7c. Approximate Steady Stale Treatment 

In order to test whether or not the role of optical modes can explain 
Ryder's data, we shall use a very crude method. We shall assume that the 
electrons all have the same energy and shall calculate their mobility on the 
basis of the mean free time at that energy; from this we calculate the power 
input. We shall also calculate the power loss in the same way. It is obvious 
that this treatment is a very poor approximation to the actual situation. 
An electron which loses energy to the optical modes will, under most cir- 
cumstances, have only a small fraction of its energy left afterwards; thus to 
assume a monoenergetic distribution is unrealistic. However, the treatment 
does bring into the analytic expressions the principal mechanisms and, as 
we shall show, appears to account for the main experimental features. 

The collision frequency or relaxation time for transitions involving the 
optical modes is given in (A7.19) and the energy loss in (A7.23). We must 
introduce corresponding expressions for the effect of the acoustical modes. 
Since the single energy distribution is to be used over the entire range of 
electric fields, we must introduce some approximations like those discussed 
in connection with (A6.25) in order to make it converge on the correct be- 
havior at E = 0. The particular choice selected is a compromise between the 
energy loss formulae for the Maxwellian and single energy distributions: 

{dZ/dl) ocous.phonons = (4 wcVOU - (vJv T y). (A7.24) 

A simplified expression is also used for the mobility: 

H = qt/mvi = noV T /v. (A7.25) 

The relationship between no and I given by this differs by 25 per cent from 
the correct relationship (A6.7); since no is an adjustable parameter in the 
comparison between theory and experiment, (A7.25) does not introduce any 
error at low fields. Equations (A7.24) and (A7.25) cause (i to converge on no 
and 81 on kT as E approaches zero. (It is probable that a slightly better fit 
to the data would be obtained by using the procedure described with equa- 
tion (A6.25); the calculations based on (A7.24) and (A7.25) were made be- 
fore the (A6.25) procedure was worked out, however, and it was not consid- 
ered worth while to rework them for this article.) 
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Rewriting the equation for mobility in terms of the collision frequencies 
1/ T = vi/{ and l/r op , the power input from the electric field is 

WdQaM = '/Mfi 2 = cfE'-fm [(1/r) + (l/r op )] 
= quotf/ivi/vr) [1 + «/4p)(»s/»i)] 
where the 7» 2 term is omitted if v t < V,. The power delivered by phonons is 
(d&/dt) phononB = 4m<* (vy/t) [1 - (vi/Vt) 2 ] ~ W4 P 

= (4^//io)(«iAr) (A7.27) 

X [1 - (vi/v T Y - (/*v/W)(/// 0P )(i>8/*i)]. 

The two coefficients of fa/«i) are both larger than unity according to the 
analysis given above. We shall introduce the symbols A and B for them: 
Accordingly 

A = l/l op = hv/kT, (A7.28) 

the last equality following from (A7.19), and 

B = hu/4 mc\ (A7.29) 

If we take hv = k 520°K, m = the electron mass and c = 5 X 10 5 cm/sec, 
we find 

B = 87. (A7.30) 

As discussed in the text, the losses appear to be larger than can be ac- 
counted for by these values of m and c. The critical drift velocity used in 
the fit of Fig. 4 was 2.6 X 10 6 cm/sec and this corresponds to a value of c 
of 

c = v /\.S\ = 1.72 X 10 6 cm/sec (A7.31) 

according to the exact treatment based on the sphere model. (As stated in 
the text this means an effectiveness of energy interchange about (1.72 X 
10 6 /5 X 10 5 ) 2 = 10 times larger than the simple theory.) 
Our simplified energy loss equation (A7.27) leads to 

v c = 2c (A7.32) 

so that we shall take 

c = 1.3 X 10 8 (A7.33) 
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in this section so as to agree with the critical velocities observed in Fig. 2. 
This leads to a value for B of 

B = k 520°K/4m (1.3 X 10 6 ) 2 

(A7.34) 
= 12.8 

For A we shall take 

A - 520/r (A7.35) 

The only other adjustable parameter is mo. For this we shall use the value 
based on Haynes' drift mobility and acoustical scattering. 

no(T) = mo (298°iO(298°A7r) 3 / 2 

(A7.36) 
= 3600 (298 K/T) 3 ' 2 

This value automatically fits the room temperature data in the Ohm's law 
range. The T~ zl - dependence then extrapolates it to the other ranges. 
The steady state condition may then be written in the form 



where 



a?(l + Ay)(ABy + Ax? - 1) - z 2 (A7.37) 

v, = (2hv/myi* (A7.38) 

x = vi/v„ y = vt/vi = (1 - .x- 2 ) 1/2 (A7.39) 

A = hv/kT = (vp/vtY (A7.40) 

B - 12.8 (A7.41) 

z = /i„ {T)E/2cA"\ ■ (A7.42) 

This form lends itself to calculation of s as a function of .r. The drift velocity 
is then found to be given by 

u = Vd/2c = s/.v(l + Ay) 

(A7.43) 
- [(ABy + Ax* -1)/(1 + Ay)] 1 ' 2 

If A » 1, there are three distinct ranges of behavior for u versus z: 

Range (/) u = z/A 1 ' 3 

For 2 — > 0, x 2 —* l/A, v = and consequently, 

u = zA 1 * 2 = v d /2c = mo E/ 2c (A7.44) 
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so that the low field relationship 

v d = noE (A7.45) 

is correctly given. 

Range II, As?s> 1 and x < 1 

In this range the electrons are at high temperature but not high enough 
to excite the optical modes. For it 

z 2 = Ax*, x = J»/AV* (A7.46) 

u = z^-A 1 '* or v d = (2cmo£) 1/2 . (A7.47) 

This corresponds to the square root range with a critical field of E e = 
2c/Vo and v e = 2c. 

Range III, x > 1 

When x is greater than unity, the optical modes enter the picture. For 
the three cases considered the values of A and AB are: 

77% A = 6.75, AB = 87, 

193°K, A = 2.69, AB = 34.5, (A7.48) 

296% A = 1.74, AB = 22.3. 

The large value of AB means that as soon as y is appreciably greater than 
zero, say 0.5 corresponding to x = 1.15, energy losses to optical modes 
dominate. As y approaches unity, the value of u is approximately 

u^AB/a+Aw (A749) 

= (hu/4mc-y' 2 /(l + A- i y>°- 

leading to 

v d (A, B) = (hv/mV'W + A-*)* 

= p,/[2(l + X- l )] l/2 . 

For the values of A and B given above, the ranges are not completely sepa- 
rated. In Fig. A2 we show the theoretical curves used in Fig. 5, together 
with the limiting lines just discussed. 

For the middle or 193°K curve, we also show the fit that would be ob- 
tained if c = 5 X 10 5 cm/sec, corresponding to B = 87 as for (A7.30). It is 
seen that this deviates much more from the data than does the curve based 
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on c = 1.3 X 10° corresponding to B = 12.8. The deviation between theory 
and experiment would be still worse at 77°K, for which temperature the 
curve of Figure hi fits the data well, as is shown in Figure 5. 




Fig. A2— The theoretical curves and their limiting forms. Some of the data from Fig. 
2 are repeated here and some additional data from the orisinal publication (Ryder and 
Shockley loc. cit.) are shown by crosses. A scale of values of x and of approximate "tem- 
peratures" is also shown. The dashed curve for T = J9.V7C is drawn for the case of the 
simple theory with c = 5 X 10 s cm/sec; the change at 77°K would be even more marked. 

Above range III, the Ax 1 term makes an appreciable contribution. 
When Ax- becomes large compared to B, the approximation of taking the 
acoustical modes to be fully excited becomes questionable. This effect may 
be estimated by comparing kT and the energy in an acoustical transition. 
The ratio is approximately 



PlC 

kT 


mxi\c hv 2c 
" If kf ' IV 




520 2-1.3-10 6 
= — • X 



(A7.51) 



300 1.26-10 7 



- x — x/3. 
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On Fig. A2 we show the values of x. Equ. (A7.51) leads to values of x < 3 
for 29S°K and x < 2 for 193°K. Although the approximation of hv (acousti- 
cal) < kT breaks down, extrapolating the acoustical scattering into the 
higher range involves some compensating effects. 

The effective "temperature" T e of the electrons may be taken on the basis 
of this approximate treatment to be proportional to v\. In order to make 
T e become equal to T for zero field, we define T e by the equation 

t. - r„!M = * V* (A7 52) 

= 520 x 2 . 

Some temperatures deduced from this equation are also shown on Fig. A2. 
Some of the data of Fig. 2 are also repeated in Fig. A2. In addition some 
earlier data 30 are also shown. These data extend to a somewhat higher 
range and appear to show the upward tendency predicted by the theory. 
The scale of temperatures indicates that for the extreme conditions experi- 
enced electron "temperatures" of about 4000°K have been produced. 

30 E. J. Ryder and W. Shockley, Phys. Rev. 81, 139 (19S0). 



